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Elliptic flow from nuclear collisions is a hadronic observ- 
able sensitive to the early stages of system evolution. We re- 
port first results on elliptic flow of charged particles at midra- 
pidity in Au-|-Au collisions at a/s^^n" ~ 1^0 GeV using the 
STAR TPC at RHIC. The elliptic flow signal, V2, averaged 
over transverse momentum, reaches values of about 6% for 
relatively peripheral collisions and decreases for the more cen- 
tral collisions. This can be interpreted as the observation of 
a higher degree of thermalization than at lower collision ener- 
gies. Pseudorapidity and transverse momentum dependence 
of elliptic flow are also presented. 

The goal of the ultrarelativistic nuclear collision pro- 
gram is the creation of a system of deconfined quarks 
and gluons [Q. If this system is created, its evolution 
should be governed by the physics of deconfined mat- 
ter. The elliptic flow observable, which is sensitive to the 
early evolution of the system, is the anisotropic emis- 
sion of particles "in" or "out" of the reaction plane de- 
fined for non-central collisions by the beam direction (z- 
axis) and the impact parameter direction (x-axis). El- 
liptic flow is usually characterized in terms of particle 
momenta by V2 — {{Px ^ pI) / {Px + Py)) ^ the second har- 
monic Fourier coefficient in the azimuthal distribution of 
particles with respect to the reaction plane [pl|| . Elliptic 
flow has its origin in the spatial anisotropy of the system 
when it is created in a non-central collision, and in parti- 
cle rescatterings in the evolving system which convert 
the spatial anisotropy to momentum anisotropy. The 
spatial anisotropy in general decreases with system ex- 
pansion, thus quenching this effect and making elliptic 
flow particularly sensitive to the early stages of the sys- 
tem evolution 0]. Being dependent on rescattering, el- 
liptic flow is sensitive to the degree of thermalization of 
the system g|[ at this early time. Hydrodynamic mod- 
els, which are based on the assumption of complete local 
thermalization, predict the strongest signals ||6|-y|. 

Elliptic flow in ultrarelativistic nuclear collisions was 
first discussed in Ref. |7| and has been studied inten- 
sively in recent years at AGS 0,|ll| and SPS |l|-|l|] 
energies. These studies have found that elliptic flow at 
high energies is "in-plane", V2 > 0, as expected from 
most models, and the pion elliptic flow for relatively pe- 
ripheral collisions increases with beam energy |1^] from 
about 2% at the top AGS energy jlO| to about 3.5% at 
the SPS [O. From transport cascade models for the full 
RHIC energy, a peak elliptic flow value of 1.5% is pre- 
dicted by UrQMD JTet c alculations and 2.5% by RQMD 
v2.4 n^ calculations ]18|| . Hydrodynamic models predict 
V2 as high as 10% [bI|S]. Details of the V2 dependence on 
beam energy and centrality are thought to be sensitive 
to the phase transition between confined and deconfined 
matter 1|Jfl,H,H,n9| (see also |ia] and references therein). 

We report here the first results on elliptic flow in 
Au-|-Au collisions at the Relativistic Heavy Ion Collider 
(RHIC) at ^/s^ = 130 GeV. The Solenoidal Tracker 



At RHIC (STAR) [|o| consists of several detector sub- 
systems in a large solenoidal magnet. For first year 
data taking, the setup consists of the Time Projection 
Chamber (TPC) which covers the pseudorapidity range 
\ri\ < 1.8 for collisions in the center of the TPC, and 
has complete azimuthal coverage, which is desirable for 
the study of azimuthal correlations. In the first year, 
the TPC is operated with a 0.25 Tesla field, allowing 
tracking of particles with pt > 75 MeV/c. Two Zero De- 
gree Calorimeters p^ which measure fragmentation neu- 
trons are used in coincidence for the trigger. The TPC 
is surrounded by a scintillator barrel which measures the 
charged particle multiplicity, and is used in studies of the 
trigger performance and vertex reconstruction efficiency. 
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FIG. 1. The primary track multiplicity distribution as a 
function of the number of tracks normalized by the maximum 
observed number of tracks. The eight centrality regions used 
in this analysis are shown. The integral under the curve is 1.0 
and the cumulative fraction corresponding to the lower edge 
of each centrality bin is also indicated in percent. 

The relative multiplicity distribution for events with a 
reconstructed primary vertex is shown in Fig. |l|. An anal- 
ysis of the trigger performance and vertex reconstruction 
efficiency, together with comparisons with Hijing p2| , 
show that the events in Fig. yj are hadronic Au+Au in- 
teractions corresponding to about 90% of the geometric 
cross section, the losses being due to vertex reconstruc- 
tion inefficiency for low-multiplicity events. This vertex 
finding inefficiency is not included in the normalization 
of Fig. y. The multiplicity is the number of primary 
tracks which pass within 3 cm of the vertex and have 
1 77 1 < 0.75. The distribution shown is not corrected for 
tracking efficiency; it is used in this analysis only to es- 
timate centralities. 

For this analysis, 22 k events were selected with a pri- 
mary vertex position within 75 cm longitudinally of the 
TPC center and within 1 cm radially of the beam line. 
Tracks were selected with 0.1 < pt < 2.0 GeV/c in order 



to have a tracking efRciency constant to within ± 10%. 
They also passed within 1 cm of the primary vertex, had 
at least 15 space points, and \r]\ < 1.3. For the determi- 
nation of the event plane we required \ri\ < 1.0. Also, the 
ratio of the number of space points to the expected max- 
imum number of space points for that particular track 
was required to be greater than 0.52, largely suppress- 
ing split tracks from being counted twice. However, the 
analysis results are not sensitive to these cuts. 
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FIG. 2. Correlation between the event plane angles deter- 
mined for two independent subevents. The upper set of data 
is for the second harmonic and the lower set is (cos(^" — '$'')) 
for the first harmonic. The points are positioned at the values 
of mean iich/nmax corresponding to each of the centrality bins 
in Fighl. The horizontal bars show the widths of the bins. 

The analysis method Ppl involves the calculation of 
the event plane angle, which is an experimental estimator 
of the real reaction plane angle. The second harmonic 
event plane angles, ^2, are calculated for two subevents, 
which are groups of independent particles from the same 
event. In order to see whether these planes are correlated, 
the mean cosine of the difference in their event plane 
angles is calculated. Although the STAR detector has 
good azimuthal symmetry, small acceptance effects in the 
calculation of the event plane angle were removed by the 
methods of shifting or weighting 0] . This correction, by 
either method, is negligible for the second harmonic. 

The subevents have been chosen in three different 
ways: 1) Assigning particles with pseudorapidity 0.05 < 
77 < 1 to one subevent and particles with — 1 < ry < —0.05 
to the other subevent. The "gap" between the two re- 
gions ensures that short range correlations, such as Bose- 
Einstein correlations or Coulomb final state interactions, 
contribute negligibly to the observed correlation. 2) Di- 
viding all particles randomly into two subevents. 3) As- 
signing positive particles to one subevent and negative 
particles to the other. Fig. g shows the results for corre- 
lation of the event planes of subevents assigned by the 
pseudorapidity method. The two other methods give 



similar results. Non-fiow effects (not correlated with the 
reaction plane) would contribute differently for these dif- 
ferent subevent choices. The shape of the centrality de- 
pendence of the second-harmonic signal is characteristic 
of anisotropic fiow and quite different from possible non- 
flow sources. 

Most commonly discussed non-flow sources of az- 
imuthal correlations are: 1) Momentum conservation, 
which can affect directed flow when each subevent is 
not symmetric about mid-rapidity, does not affect elliptic 
flow measurements. 2) Coulomb and Bose-Einstein cor- 
relations p3| , which are eliminated by the construction 
of the subevents in Fig |[ 3) Resonance decay Q , whose 
effect on the subevent correlation would be independent 
of centrality, unlike what is observed for the second har- 
monic. 4) Jets, when calculated using Hijing ||2^ for the 
cuts used in the current analysis, do not contribute be- 
yond the systematic errors for V2 quoted below. Also, if 
jets or resonances contribute to V2 they would contribute 
to the directed ffow measurements in comparable amount 
|25| . The first harmonic correlation, which is shown in 
Fig g, and the higher harmonics, are about an order of 
magnitude weaker than the elliptic fiow at mid-centrality. 
This sets an upper limit for the contribution of all non- 
flow effects to the elliptic flow and is the basis of the 
estimate below of the systematic errors. 

The analysis method involves correlation of the az- 
imuthal angle, 0, of each particle with an event plane 
angle, ^, and then averaging over all events. In this pa- 
per we have used three particle correlation methods: 1) 
Correlating the particles from one hemisphere with the 
event plane of the subevent in the other hemisphere. 2) 
Correlating each particle with the event plane of all the 
other particles. 3) Correlating particles of one charge 
sign with the event plane of the opposite charge sign. 
The observed elliptic flow comes from the second har- 
monic Fourier coefficient of the particle azimuthal distri- 
bution with respect to the event plane, which is simply 
{cos[2{(j) — ^2)])- The elliptic flow relative to the real 
reaction plane, '^r, the plane defined by the impact pa- 
rameter and the beam, can be evaluated by dividing the 
observed signal by the resolution, (cos[2(^ — ^a)]), of the 
event plane. The resolutions calculated from the correla- 
tion of subevent planes were somewhat different for the 
different subevent selections, but the resultant V2 values 
were the same within statistical errors. The resolutions 
for the full events reach 0.7 for the centrality at the peak 
in Fig. |, while in NA49 |l|,|l| at the SPS they only 
reached 0.4. A resolution of the event plane angle of 0.7 
is sufficiently close to the ideal value of 1.0, to correlate 
other quantities, such as two-particle correlation mea- 
surements (HBT) with the event plane. Since we do not 
measure the correlation with the first harmonic plane, we 
cannot determine the sign of V2 ■ 

Our analysis procedures have been tested with simu- 
lated data pa to which a known amount of fiow has been 



added. The simulated data were filtered by a GEANT 
model of STAR and reconstructed in a way similar to that 
used for the data. For 2% and 10% elliptic flow added to 
the simulations, the flow extracted was (2.0 ± 0.1)% and 
(9.7 ± 0.2)%, respectively. 
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FIG. 3. Elliptic flow (solid points) as a function of central- 
ity defined as rich/nmax- The open rectangles show a range of 
values expected for V2 in the hydrodynamic limit, scaled from 
e, the initial space eccentricity of the overlap region. The 
lower edges correspond to e multiplied by 0.19 and the upper 
edges to e multiplied by 0.25. 

Fig. H shows V2 as a function of centrality of the colli- 
sion. Although this figure was made with the subevents 
chosen as in Fig. 0, the same results within errors were 
obtained with the other correlation methods. Restricting 
the primary vertex z position to reduce TPC acceptance 
edge effects also made no difference. From the results of 
the study of non-flow contributions by different subevent 
selections and the maximum magnitudes of the first and 
higher-order harmonics, we estimate a systematic error 
for V2 of about 0.005, with somewhat smaller uncertainty 
for the mid-centralities where the resolution of the event 
plane is high. The systematic errors are not included in 
the figures. 

In the hydrodynamic limit, elliptic flow is approxi- 
mately proportional to the initial space anisotropy, e, 

1 27]]. The transformation to 
! was done using a Hijing |2|] 



which is calculated in Rcf. 
the multiplicity axis in Fig. I 
simulation, taking into account the above mentioned 
vertex- finding incfhcicncy for low multiplicity events. In 
comparing the flow results to e, no unusual structure is 
evident which could be attributed to the crossing of a 
phase transition while varying centrality H,M. The e 
values in Fig. ^ are scaled to show the range of hydrody- 
namic predictions ^^ for V2/e from 0.19 to 0.25. The 
data values for the lower multiplicities could indicate in- 
complete thermalization during the early time when el- 
liptic flow is generated |^,^. On the other hand, for 
the most central collisions, comparison of the data with 
hydrodynamic calculations suggest that early-time ther- 



malization may be complete. The V2 values peak at more 
peripheral collisions than RQMD predictions ||T^ , but in 
qualitative agreement with hydrodynamic models [Q. 
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FIG. 4. Elliptic fiow as a function of transverse momen- 
tum for minimum bias events. 

The differential anisotropic flow is a function of 77 and 
Pt- For the integrated results presented here, all v val- 
ues should first be calculated as a function of rj and pt, 
and then averaged over either or both variables using the 
double differential cross sections as weights. Since we do 
not yet know the cross sections, we have averaged us- 
ing the observed yields. Fig. H shows V2 as a function of 
Pt for a minimum bias trigger. The ry dependence (not 
shown), which is averaged over pt from 0.1 to 2.0 GeV/c, 
is constant at a value of (4.5 ± 0.5)% for \r]\ < 1.3. We 
have assumed that the efficiency (yield/cross section) is 
constant in the pt range where the yield is large. This is 
borne out by studies of the effects of different track qual- 
ity cuts on the observed pt spectra. For the pt depen- 
dence the data are not very sensitive to the assumption 
of constant efficiency as a function of rj because V2 ap- 
pears to be independent of rj in the range used, \r]\ < 1.3. 
Mathematically the V2 value at pt = 0, as well as its 
first derivative, must be zero, but it is interesting that V2 
appears to rise almost linearly with pt starting from rela- 
tively low values of pt ■ This is consistent with a stronger 
"in-plane" hydrodynamic expansion of the system than 
the average radial expansion. Note that the results shown 
in Fig. @ were obtained by taking the average over both 77 
and Pt, weighted by the yield. Although Fig. ^ is for ap- 
proximately minimum bias data [^ the general shapes 
are the same for data selected on centrality, except that 
the slopes of the pt curves depend on centrality. Fig. ^ 
was made using pseudorapidity subevents, although the 
same results within errors were obtained using the other 
two methods. 

We conclude that elliptic flow at RHIC rises up to 
about 6% for the most peripheral collisions, a value which 



is more than 50% larger than at the SPS |l^, indicating 
stronger early-time thermalization at this RHIC energy. 
Elliptic flow appears to be independent of pseudorapidity 
in the region \r]\ ^ 1.3. Its pt dependence is almost linear 
in the region 0.1 < pt < 2 GeV/c. Comparing to esti- 
mates @ based on transport cascade models, we find 
that elliptic flow is underpredicted by RQMD by a fac- 
tor of more than 2. Hydrodynamic calculations |^,|[ for 
RHIC energies overpredict elliptic flow by about 20-50%. 
This is just the reverse of the situation at the SPS where 
RQMD gave a reasonable description of the data and hy- 
drodynamic calculations were more than a factor of two 
too high [0|. Also in contrast to lower collision energies, 
the observed shape of the centrality dependence of the 
elliptic flow is similar to hydrodynamic calculations and 
thus consistent with signiflcant thermalization. The val- 
ues for elliptic flow compared to hydrodynamic models 
indicate that early-time thermalization is somewhat in- 
complete for peripheral collisions but perhaps complete 
for the more central collisions. 
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